Hadronic Effective Field Theory Applied to A-Hypernuclei 

Jeff McIntire 

Department of Physics, College of William and Mary, Williamsburg, VA 23187 



In the present work, the approach of Furnstahl, Serot, and Tang (FST) 
is extended to the region of nonzero strangeness in application to single- 
particle states in single A-hypernuclei. To include A's, an additional con- 
tribution to their effective lagrangian is systematically constructed within 
the framework of FST. The relativistic Hartree (Kohn-Sham) equations are 
solved numerically, and least-square fits to a series of experimental levels are 
performed at various levels of truncation in the extended lagrangian. The 
ground-state properties of any A-hypernuclei are then predicted. In addi- 
tion, ground-state A-particle-nuclcon-holc splittings are calculated where 
appropriate, and the approach is calibrated against a calculation of the 
Si/2-doublet splitting in the nucleus IfPn- 

PACS numbers: 21.80.-|-a 



1. Introduction 

Effective field theories have been developed in recent years to solve the 
nuclear many-body problem. In the present work, we consider one of these 
theories, proposed by Furnstahl, Serot, and Tang (FST) ^ [2], and ex- 
tend their methodology to the case of single A-hypernuclei. Specifically, the 
phenomena of interest here are ground-state (GS) binding energies, densi- 
ties, single-particle spectra, and particle- hole splittings of select single A- 
hypernuclei. 

FST develop a self-consistent framework for constructing an effective 
lagrangian that incorporates the principles of both quantum mechanics and 
special relativity, the underlying symmetries of QCD, and the nonlinear real- 
ization of spontaneously broken chiral symmetry 1 . As this is a low-energy 
theory, the appropriate low lying hadrons are used as degrees of freedom. In 
order to make any meaningful calculation, the lagrangian, which in principle 
contains an infinite number of terms, must be truncated in some way. Naive 
dimensional analysis (NDA) and relativistic mean field theory (RMFT) 
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[21 E] provide a formalism in which higher order terms are, in general, suc- 
cessively smaller; this allows for a systematic expansion, and consequently 
a meaningful truncation, in the effective lagrangian. Here FST utilize rela- 
tivistic Hartree theory to reduce the many-body equations to single-particle 
equations. The free parameters in the effective lagrangian are fixed via 
least-squares fits to experimental data from ordinary nuclei along the valley 
of stability. These fits are conducted at various levels of truncation in the 
underlying lagrangian Once the values of these parameters are known, 
this lagrangian can be used to predict other properties of ordinary nuclei. 
One example which demonstrates the predictive power of this method is its 
application to the study of nuclei far from stability [H| |7j . 

Density functional theory (DFT) is a theoretical framework which al- 
lows one to calculate the GS properties of many-body systems without car- 
rying around all the baggage contained in the many-particle wave functions 
[S]. Two points are of interest here. First, if the expectation value of 
the hamiltonian is considered as a functional of the density, the exact GS 
density can be determined by minimizing the energy functional. Second, 
one only needs to solve a series of self-consistent, single-particle equations 
with classical fields, instead of many-body equations with quantum fields 
[S]. In other words, Kohn-Sham theory is formally equivalent to relativistic 
Hartree theory. Consequently, the problem is now reduced to determining 
the correct form of the energy functional, which follows from the appropri- 
ate lagrangian. The full interacting lagrangian of FST gives an appropriate 
energy functional and, as a result, DFT provides an underlying theoretical 
justification for this approach. 

Hadronic effective lagrangians using MFT have been developed in the 
literature to describe hypernuclei. Early models containing only the lowest 
order terms required much weaker meson couplings to the A than to the nu- 
cleons to achieve success |1U1 lllj , particularly in the weak spin-orbit interac- 
tion. Later, it was suggested that large meson couplings to the A consistent 
with SU(3) were possible if the lagrangian was extended to include tensor 
couplings dUISlITllITSlIinilllllTHI. It turns out the spin-orbit splitting is 
very sensitive to the size of the tensor coupling to the vector field. The ap- 
proach of FST has also been applied to strange hadronic matter ^H]- More 
recently, effective theories consistent with SU(3)l <8) SU(3)r, have been con- 
structed [201 [^ [221 • Another model of interest uses strangeness changing 
response functions to calculate the spectra of yO and yCa;^ the resulting 
GS particle-hole splittings are small [23 ■ Other studies include models that 
couple the mesons self-consistently to the quarks within the baryons [241 125j 
and a density dependent relativistic hadronic field theory [26j . 



^ Here Y denotes a hyperon. 
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The following studies have attempted to fit potentials to the hyperon- 
nucleon interaction. Experimental data has been analyzed to obtain a non- 
local and density-dependent A-nucleus potential [^Tj |2H]- Global optical 
potentials for A scattering off nuclei were developed j^HI- The hypernuclear 
mass dependence of the binding energies is reproduced by a A moving in a 
Woods-Saxon potential [S^. The Nijmegen group has developed Y-N po- 
tentials based on the assumption of SU(3) symmetry [^OHISS]; this fixes 
the baryon-meson coupling constants from N-N scattering fits. Similarly, 
potentials were constructed by the Julich group assuming SU(6) symmetry 
j34j . Calculations of hypernuclei using these Nijmegen or Julich poten- 
tials include [SHI OHl IHZl IHHl IHSl 1^- Comparable G-matrix calculations 
with a SU(6) quark-model baryon-baryon interaction |41| and Skyrme-like 
hyperon-nucleon potentials have also been investigated. Other recent 
approaches include using the Fermi hypernetted chain method [431 144j and 
using a quark model with one boson exchange potentials |45j . 

Many of these studies achieve a good deal of success. However, the 
framework of FST is more comprehensive than these approaches as it in- 
corporates directly into a hadronic effective field theory all of the following: 
special relativity, quantum mechanics, the underlying symmetry structure 
of QCD, and the nonlinear realization of spontaneously broken chiral sym- 
metry. Furthermore, this methodology is theoretically justified by DFT. 
Once all the parameters are fixed, their lagrangian predicts the GS proper- 
ties of any ordinary nucleus. This approach has had great success 1, 61 [7|. 
Therefore, it is of interest to extend this methodology, with all of its intrinsic 
strengths, to the strangeness sector, as is done here. 

In the present work, the approach developed by FST is expanded to 
the region of the strangeness sector that corresponds to A-hypernuclei with 
S = — 1 and T = 0. To this end, we include a single, isoscalar A field in the 
theory.^ Now, a A-lagrangian is constructed as an additional contribution 
to the full interacting effective lagrangian of FST, consistent with their 
methodology. Since the A is an isoscalar, it does not couple to either a 
single Yukawa pion or the rho meson. Furthermore, we confine our theory 
to the mesons already included;^ thus, the meson lagrangian, which in this 
approach contains the majority of the complexity, is unaltered. It has been 
proposed that a tensor coupling to the vector field be included to reproduce 
the correct experimental spin-orbit splitting of the p-states in A-hypernuclei 



^ The E is not explicitly included in the present calculation. An idea of the possible 
impact of A — E mixing can be taken from |46) . It should be mentioned that if one 
views the scalar meson as a two-pion resonance, then the E enters implicitly as an 
intermediate state in our formalism. 

The kaon is not included as a degree of freedom in this work. The reason is that, as 
with the pion, the kaon has no mean field and does not effect the RMFT calculations. 
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jl2|ll,Sj. As it turns out, such a term is a natural extension of our lagrangian 
in this framework. Additional higher order terms are also included to better 
approximate the exact energy functional. 

Following the methodology of FST, our A-lagrangian contains a number 
of free parameters. The constants in both the nucleon and meson sectors 
are taken from a FST parameter set corresponding to their full lagrangian. 
As before, the remaining unconstrained parameters are fixed here via least- 
squares fits to a series of experimental data [13 UHl EHl EOl lEIl IHlj' The 
10 pieces of data used include six GS binding energies, three s-p shell ex- 
citations for the A, and the spin-orbit splitting of the p-states in j^C. The 
fits are conducted at four different levels of truncation in the A-lagrangian. 
Once these parameters are fixed, this lagrangian can be used to predict 
other properties of single A-hypernuclei. 

One other property that is of interest to calculate here is what we refer 
to as Si/2-splittings. These are GS particle-hole splittings of select single 
A-hypernuclei, such as j^O, which have a A in the GS and a hole in the 
last filled nucleon (proton or neutron) shell. For these systems, the angular 
momenta of the A and the nucleon hole couple to form a doublet. The 
size of these splittings is determined by the difference of two particle-hole 
matrix elements . The effective particle- hole interaction utilized here 
follows directly from the effective theory of the preceding discussion. This 
interaction, to lowest order, is just simple scalar and vector meson exchange 
j54j.'^ A simple Yukawa spatial dependence is obtained when retardation is 
neglected in the meson propagators. With this exception, the full Lorentz 
structure is maintained jSl]. For the A-N case, there is no isovector compo- 
nent to the effective interaction or exchange contribution in the two-body 
matrix elements. Through angular momentum relations j57j and some al- 
gebra, the matrix elements are reduced to radial Slater integrals. Using the 
Hartree wave functions from the A single-particle calculations to evaluate 
the integrals, these matrix elements, and consequently the Si/2-splitting, 
can now be fully determined. Once the parameters in the A-lagrangian are 
known, the effective particle-hole interaction is completely specified in this 
approach. In the case of Si/2"Splittings in A-hypernuclei, only the spatial 
part of the vector exchange contributes to the splitting. Predictions are 
given for the GS doublet splittings of every one of the A-hypernuclei con- 



* The retention of higher diagrams in the effective interaction, particularly those in- 
cluding the tensor coupling to the A, is left for future work. Also, it is worth noting 
that while the kaon makes no contribution at the mean field level, kaon exchange may 
play a role in the effective interaction. Some idea of the relative contribution of kaon 
exchange can be obtained from the Nijmegen potentials '32' SKIIS^ . An investigation 
of the effect of kaon exchange on the Si/2-splittings in effective field theory is also left 
to future work. 
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sidered here; all of the doublets used in the fitting procedure lie within the 
current experimental error bars on the GS energies. An upcoming high res- 
olution experiment at Jefferson Lab will measure the Si/2-splittings in j^^B 
and j^^N [HHllEn]- The present calculations provide theoretical predictions 
for these quantities!^ Non-relativistic calculations of similar particle-hole 
splittings have been carried out |6n] . 

The need for isovector interactions and exchange contributions make 
calculations of similar splittings in ordinary nuclei far more complicated 
j54j . As an example of a comparable system in an ordinary nucleus, and 
to at least partially calibrate the present approach, the calculation of the 
Si/2-splitting in ffPi? is included here. Comparable systems for ordinary 
nuclei have also been examined |61j . 

In section |^ we review the methodology of FST and in section |21 we 
describe the development of our A-lagrangian. The framework for calculat- 
ing the particle-hole splittings is discussed in section |3 The results of the 
parameter fits, single-particle calculations, and Si/2-splittings are given in 
section 13 



2. Methodology of FST 

In this section we review the methodology of FST. They approach the 
nuclear many-body problem by constructing an effective field theory that 
retains the underlying symmetries of QCD as well as the principles of both 
special relativity and quantum mechanics At low-energy, hadrons are 
the desired degrees of freedom and the ones which FST use to construct 
an effective lagrangian. The nonlinear realization of spontaneously broken 
chiral symmetry is illustrated through a system of pions, nucleons, and rho 
mesons. They incorporate Goldstone pions through the field 

U(x^) = C(x^)ie(x;.) = e^-(-M)//.ie*-(-M)//. (i) 

where the pion field, 7r(x^) = • tt, appears to all orders, r is a Pauli 
matrix, and f^^ is the pion-decay constant PP. An isodoublet nucleon field 
is included, represented by 

^<^^) - ( ) 

The upper (lower) component corresponds to the proton (neutron). To 
account for the symmetry energy in nuclear matter, an isovector-vector rho 
meson, Pui'^/i) = ' P, is also included. 

The following boson fields are also incorporated into this framework, the 
first two of which are isoscalar chiral singlets. A scalar field, (j), is included 
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to simulate the medium-range nuclear attraction. Next, they incorporate a 
vector meson, V^, to reproduce the short-range nuclear repulsion. Lastly, 
a photon field, A^, is added to take into consideration the electromagnetic 
structure of nuclei. 

As all possible combinations of the fields, consistent with this frame- 
work, are included, this lagrangian contains an infinite number of terms. 
To conduct any meaningful calculation, this lagrangian needs to be trun- 
cated at some level. FST utilize both NDA and RMFT to accomplish this. 
NDA is a framework which identifies all the dimensional factors of a given 
term. Once these dimensional factors, and some appropriate counting fac- 
tors, are extracted from a term, the remaining dimensionless constant is 
of 0(1) This assumption is known as "naturalness." RMFT states 

that when the baryon density becomes appropriately large, the sources and 
meson fields can be replaced by their expectation values; here, the expec- 
tation values of the meson fields are just their classical fields 0- Then we 
notice that while the meson mean fields are large, the ratios of these fields 
to the chiral symmetry breaking scale, M, are small. Furthermore, the size 
of derivatives is related to k-p, which is also small compared to M. These 
effects are shown by 

^ W 1 kp 1 

M ' M 3' M 4 ^ ^ 

where the scaled meson mean fields are defined as 



$(x) = gs(/)o; W(x) = gvVo; R(x) = gpbo; A(x) = eAo (4) 
The ordering principle developed by FST is 



n 



+ b + d 



(5) 



where for a given term v is the order, n is the number of fermion fields, h 
is the number of non-Goldstone bosons, and d is the number of derivatives. 
Now a controlled expansion is performed in which higher order terms are, 
in general, progressively smaller. 

Using this ordering principle, they construct an effective lagrangian in 
two parts jJI 

>Cfst(x^) = >Cn(x^) + >Cm(x^) (6) 
The fermion part to order = 4 is given by 

d 



dx, 
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|eA^ (1 + Ts) 



+ (M-gs</.)|N + ^NcT^,p^,N 



1G — U 

+ ^N>(/^s + /?vr3)N— F,. (7) 

where A = iAp(l + T3) + ^Xn{l - T3) and Ap = 1.793 (An = -1.913) is the 
anomalous magnetic moment of the proton (neutron). Note that for the 
purposes of this work, the conventions of fS^ are used. Here we have defined 

■W/x!/, Pfiu, and F^i, are similarly defined for v^, p^, and respectively. 
Notice that the pions only couple to the fermions through the combinations 



(^) 



2 \ (?x^ dxf^ 



To lowest order, both and contain derivatives of the pion field; thus 
soft pions decouple. The meson lagrangian to order = 4 is 

- - ( 1 + I ^t^U^t^U - gpTTTT^tr [Pf^uVf^u) 

.^tr(u.Ut-2)-l(l..f .|f!)m^V.V. 
+ ^Cog^ (V^V^)2 - If^^F^, - (^1 + r?,|^) m^tr (/,^p^) 

Terms such as NNc/)^ are redundant in this formulation. This stems from 
the fact that FST employ meson field redefinitions; since the parameters 
are free, they are also just redefined. A detailed description of how this 
lagrangian was constructed is presented in pp. 
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This still constitutes a system of many-body equations with quantum 
fields. FST now employ Hartree theory and RMFT to reduce the many- 
body system to a series of single-particle equations with classical fields. This 
is equivalent to Kohn-Sham theory in DFT; therefore, DFT provides the 
theoretical justification for this methodology. The single-particle hamiltonian 
takes the form fT| 

h(x) = -ia- V + W + ^rsR + ^ (1 + rs) A + /3 (M - gs«>) - ^A/3a- VA 



1 



1 



—13a- [ /vVW + /p-r3VR ) + ^ (/3s + ^^r^) V'A 



2M2 



(12) 



Since the pion has no mean field in a spherically symmetric system, all of 
the pion couplings drop out. The Hartree wave functions are of the form 



V'a(x) = - 



fGa(r)$«in 



Ct 



(13) 



Here a = {a, m} = {nlsj,m}, is a two component spinor, and ta is 1/2 
(-1/2) for protons (neutrons). The are the spin spherical harmonics. 
Substituting this wave function into the Dirac equation, 



h(x)V'a(x) = EaVa(x) 

one arrives at the following radial Hartree equations 
d 



K 

dr r 



" d K 
_di r 

where the single-particle potentials are 



Ga{r) - [E, - Ui + U2] F,(r) - UgGJr) = 
Fa(r) + [E, - Ui - U2] Ga(r) + U3F,(r) = 



(14) 

(15) 
(16) 



Ui(r) = W(r) + t,R(r) +[ta + -] A(r) + 



1 



1 



U2(r) = M-'l'(r) 



+ 



1 gA(r; 
2M di 



(Ap + An) 



2 

aR(r 
di 

+ ia(Ap — An) 



2M2 



(/3s + 2t,/3v) V2A(r) (17) 
(18) 



(19) 



The scalar meson equation is determined by minimizing the variational 
derivative of the effective lagrangian with respect to the scalar meson field. 
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The other meson equations are constructed in a similar fashion. These 
meson equations are P 



V^^ + m^$ = gips(x) -^^'[^ + ^-, 



3! M 



si 



2M 



■M 



7^^ Ul+^?2T7 ^W^ + -i- 



r2 



2M 



2Mgv 

2ijr I ™2 TTr „2 



-V^W + m^W = gy 



2M g2 
ri2^\ ^ 2 



(20) 



Pb(x) + ^V. (pT(x)f 



+ 



02 

M 



3g7 



-V^R + m^R 



2 Bp 



«(x) + 4v.(pjwf) 



Pchg(x) 



^P^m,R 



/5chg(x) 



-V^A = eVchg(x) 



(21) 

(22) 
(23) 



The baryon sources become the densities in the meson equations and are 
given here by 



Ps(x) 


a 


(24) 


/'b(x) 


a 


(25) 


pg(x) 


= E^2G.(r)F.(r) 

a 


(26) 


/'six) 


a 


(27) 


/^J(x) 


= E^|^(2ia)2G,(r)F,(r) 


(28) 



The charge density is made up of two components 

Pchg(x) = Pd(x) + /Om(x) 



(29) 
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where the first, the direct nucleon charge density, is 

1 



Pd 



:x) = /Op(x) + ^V- (pj(x)f) + 



2M2 



/3sVVb + /?vVV3 



and the second, the vector meson contribution, is 
Pm(x) = 



— v2r+ .^^V^W 



&76p 3g^gv 

The point proton and nucleon tensor densities in Eq. H3U|) are 
^pW = ^E^^(l + 2*a)(G^(r) + F^(r)) 

a 

= ^(PB+Ps) 

/^a(x) = E^^2AG.(r)F.(r) 

a 

respectively. Finally, the energy functional is given by ^ 



(30) 



(31) 



d^xUr 



(32) 
(33) 

(34) 



where 



-^ps + -W ( PB + ^ V. p^i ]+-R{p3 + ^ V- f ) + -Ap, 



2M 



+ 



$3 



gi M 
1 



11 _L 

12 24 M 



4 M g^ 



-— ??i +172 — 

' -2 4M V M 



m- 



V 



W2 



^ (VW)^ 



4g^M 



(35) 



The radial Hartree equations and the meson equations form a system which 
is solved self-consistently until a global convergence is achieved. FST wrote 
a program to numerically solve the coupled, local, nonlinear, differential 
equations. Huertas has written an independent program which reproduces 
the results of FST [HI E] . The free parameters in this system are listed in 
Table ^ These are fit by FST to a series of experimental data along the 
valley of stability at various levels of truncation in the underlying effec- 
tive lagrangian pLj. The last three parameters are fit to the electromagnetic 
properties of the nucleon. The remaining constants are determined by mini- 
mizing a least-squares fit where 29 pieces of experimental data were used. 
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ms/M 


gs/4vr 


gv/4vr 


gp/4vr 


^1 


^^3 


G2 


0.55410 


0.83522 


1.01560 


0.75467 


0.64992 


3.2467 






/v/4 


/p/4 


m 




Co 


G2 


0.3901 


0.1734 


0.9619 


0.10975 


0.63152 


2.6416 




/3s 


/3v 


ai 








G2 


-0.09328 


-0.45964 


1.7234 


-1.5798 







Table 1. The G2 parameter set developed by FST 1 . The first 4 parameters 
correspond to v = 2^ the next 5 to = 3, the following 5 to = 4, and the last 2 
to V — 'b. 



The result of a parameter fit corresponding to their full lagrangian is shown 
in Table ^ Note that these parameters do indeed satisfy the naturalness 
assumption made earlier and as a result, higher order terms are successively 
smaller. Also, we mention that increasing the level of truncation beyond 
that of the G2 parameter set does not significantly improve the fit . Once 
the free parameters are determined, this lagrangian can he used to predict 
other properties of ordinary nuclei ^IHl CI • 



3. Application to A-hypernuclei 

We now consider an extension of this approach to the strangeness sector. 
The specific phenomena that we seek to investigate here are GS binding 
energies (i.e. chemical potentials), densities, single-particle spectra, and 
particle-hole states of single A-hypernuclei. To this end we add a single, 
isoscalar A to the theory. Note that the A is also a chiral singlet. Then, 
we construct our effective A-lagrangian as an additional contribution to the 
full v = 4 lagrangian of FST, utilizing their methodology. This lagrangian 
is of the form 



(36) 



Here we restrict ourselves to the mesons already incorporated into the theory 
by FST; therefore, the A-lagrangian is confined to the fermion sector. First, 
we consider all possible contributions up to order v = 2, consistent with 
this approach. Our effective A-lagrangian now takes the form 



C 



(2) 



-A 



7a. 



_d_ 



«gvAV^ + (Ma - gsA<; 



A 



(37) 



Notice that the coupling constants, gsA and gvA, are free parameters and 
are different from those used in the nucleon case. Single Yukawa rho and 
pion couplings to the A are absent as they do not conserve isospin. Also, no 
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electromagnetic coupling is retained to this order as Q = for the A. Four 
fermion terms are discussed in appendix 1X1 

However, this lagrangian, to order v = 2, fails to reproduce the small 
experimental spin-orbit splitting of the p-states, as in j^C |5Uj . It was 
proposed in the literature that tensor couplings of order v = She introduced 
to correct for this limitation . We add tensor couplings to the vector 

and photon fields, shown by 

^(T) ^ g^^^^^v^^^ + _^AA^^^,F^,A (38) 

The constant gxA is a free parameter. Here Aa = —0.613 is the anomalous 
magnetic moment of the A. Since we want to make a full expansion in our 
A-lagrangian to order v = 3, consistent with this approach, we must also 
include three additional terms, shown by the following 

^^^^ = ^'Wl^^^' + ^'WL^^^''^'' + i/i3^A7^A<^V^ (39) 

where ^ii, fX2, and are three more free parameters. In the nucleon case, 
the terms comparable to these last three were regrouped through redefinition 
of the meson fields. However, in the A case this is no longer possible unless 
additional mesons are added to the theory. A more complete description of 
how the terms in the A-lagrangian are chosen is contained in appendix \K\ 
Now our A-lagrangian, complete to order = 3, is 

U = C^l^+C^P+C^P (40) 

Note that our lagrangian in Eq. H36|) includes all possible terms up to = 4 
in the nucleon and meson sectors as well.^ 

In the Hartree formalism, we add a new wave function for each new 
baryon, given here for the A by 

^= 7 I -FA(r)c^.«. ) (41) 

^ It is of potential interest to consider coupling additional scalar and vector mesons, 
such as the fo and the $, to the strangeness density and conserved strangeness cur- 
rent respectively. This allows one to eliminate the terms in -C^' using the equations 
of motion and redefinitions of the new fields. However, the number of additional 
terms, and their accompanying free parameters, introduced to v = 2> make this ap- 
proach more complex than the present framework. Fortunately, the point is relatively 
unimportant for the single A-hypernuclei considered here as these new mesons are 
self-fields of the A. If they are included, they would appear only in the energy func- 
tional and have no effect on the energy eigenvalues; as the last eigenvalue in this 
approach is equivalent to the total binding energy per baryon for the GS, they have 
no efi^ect on the cases of interest here. 
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Plugging this wave function into the Dirac equation yields the following new 
pair of Hartree equations 



d K 
di r 



GA(r) - [Ea - U4 + Us] FA(r) - UgGAlr) 
FA(r) + [Ea - U4 - Us] GA(r) + UgFAlr) 







" d K 
dr r 

where the A single-particle potentials are 



U4 

Us 
Ufi 



gVA 

gv 
Ma 



W- 



gSA 



/f3 
M 



gs 

gTA dW 



+ 



2M 
Aa dA 



2M 



(42) 
(43) 

(44) 
(45) 
(46) 



2M di 2M dv 

Since all our additional terms are in the fermion lagrangian, the only change 
to the meson equations are added contributions to the source terms. The 
new contributions to the source terms arising from the A-lagrangian are 



6pi 



1 



47rr2 
1 

47rr2 
1 

47rr2 
1 

47rr2 
1 

47rr2 
1 

47rr2 



Ga 



(r) - Fi(r)) 



gs M 



Gi(r) 



+ Fi(r) 



(Gi(r) + Fi(r) 



M 



?VA 



?v 



M 



Gi(r) 



Fi(r) 



/i2 

M 



W 



2GA(r)FA(r) 



gTA 

/v 

2AAGA(r)FA(r) 



(47) 

(48) 
(49) 
(50) 



The new energy functional is identical in form to the one used by EST, 
with only one additional energy eigenvalue, Ea- The numerical solution 
to the extended set of coupled, local, nonlinear, differential equations was 
obtained by extension of a program developed by Huertas [HllZl- Here we 
use the parameter sets of EST for the nucleon and meson sectors. There 
are six new parameters in our A-lagrangian: gsAj gVA, gTA, P-i, P2, and 
//3. Least-squares fits to a series of experimentally known A single-particle 
levels are conducted at various levels of truncation in our A-lagrangian, 
while maintaining the full lagrangian of EST to order i/ = 4. Now this 
lagrangian can be used to predict other properties of single A-hypernuclei. 
One application we investigate in the next section is Si/2-splittings. 
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4. Si/2-doublets 

Consider nuclei like ]^0; the GSs of such systems are particle-hole states. 
One process by which nuclei of this type are created is the reaction (Tr"*", K~^) 
on target nuclei with closed proton and neutron shells |47| I48[ I49j . During 
the course of this reaction a neutron is converted into a A. As a result, 
a neutron hole is also created which, for the GS, inhabits the outermost 
neutron shell. The angular momentum of the A and the neutron hole couple 
to form a multiplet. Since the A occupies the lsi/2 shell in the GS, there 
are only two states in these multiplets. It is these configurations that we 
refer to as Si/2-doublets. The reaction (e, e'K"^) is another process used to 
create nuclei of this type [581 159j . This process differs in that a proton hole 
is created here and that greater resolution is possible. 

In order to calculate the splitting of these doublets, we first consider 
Dirac two-body matrix elements of the forms j54j 

((ni/iji)(n2/2j2)JM|V(ri2)|(n3/3j3)(n4Z4j4)J'M') (51) 

and 

((ni/iji)(n2/2j2)JM|V(ri2)<T(^) • a(')|(n3/3j3)(n4/4j4)J'M') (52) 

where the single-particle wave functions are specified by {n/j}, correspond- 
ing to either the upper or lower components in Eq. (|l,3j) . and V(ri2) is some 
effective interaction. Next, we expand this effective interaction in terms of 
Legendre polynomials j54j 

oo 

V(ri2) = E/fc(ri'^2)Pfc(cos^i2) 

k=0 

oo 

= E/'^(ri'^2)Cfe(l)-Cfc(2) 

A:=0 

where [HZ] 

^k,= [^kTT) ^k,iOA) (55) 
Inverting Eq. (j58|) yields the expression 

2k -\- 1 

/fc(ri,r2) = J ^d(cos0i2)Pfc(cos0i2)V(ri2) (56) 

In the case of Eq. (|52|) . the effective interaction is coupled to Pauli matrices. 
Therefore, Eq. is modified to 

Y{n,)aW . a(') = j:i--^)'^'~'Mn,r2)xt'\l) ■ xt'\2) (57) 

k\ 



(53) 
(54) 
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Here Xx^'^ are Ckq coupled to Pauli matrices, shown by 



(fc,i) 



(58) 



Now we introduce a specific type of effective interaction. The form 
we use here follows directly from the effective lagrangian in the preceding 
section and to lowest order, corresponds to simple Yukawa couplings of both 
the scalar and vector fields, given by 



V(ri2) 



(1) (2) 
74 74 



-gsgSAe ^^^'^ .,(i)./2) gvgVA e- 



47r 



ri2 



mvri2 

ri2 



(59) 



Here ri2 = |ri — r2|. This simplistic spatial dependence is possible because 
retardation in the meson propagators is neglected, or p^ = (p, P4) — > (p, 0). 
Otherwise the full Lorentz structure is maintained 54 . Couplings to the 
rho and pion fields are absent as T = for the A. In this formalism, we can 
now write 



where 



/fc(ri,r2) 



/f(ri,r2) 



(1) (2) 

71 74 



gSgSA,„, , ,N2ms. , - , ^ 
-{2k + 1) ifc(msr<)kfc(msr>, 



vr 



?VgVA.„, , , ^ 2mv . 

— {2k + 1) ifc(mvr<)kfc(mvr>; 

47r vr 



(60) 



(61) 
(62) 



where r< (r>) is the smaller (larger) of ri and r2. Here ifc(mr) and kfc(mr) 
are modified spherical Bessel functions of order k. 

The matrix elements in Eqs. (|5H) and (|52j) are actually six dimensional 
integrals. Treating the 7-matrices as 2 x 2 block matrices operating on 
the upper and lower components of the Hartree spinors, these Dirac matrix 
elements, for each term in the interaction, are actually the sum of four 
separate integrals. The scalar and vector time (/i = 4) components of the 
effective interaction take the form of Eq. H51|) : the vector spatial (/i = 1,2,3) 
components take the form of Eq. ()52() . Thankfully, angular momentum 
relations allow one to integrate out the angular dependence |57j. These 
integrals, for the scalar and vector time components, become 



(51)= 5:(12|/^(ri,r2)|34)(-l)j^+J 
fc=0 



3+J /J j2 jl 

k h u 



^jj"5mm' 



x((/4)jl||C,(l)||(/3^)j3)((/2^)j2||C,(2)||(/4^)j4) (63) 
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where i = S,V and (51) indicates the quantity in Eq. (|51j) . For the vector 
spatial components, these integrals become 

oo 

(52) = ^^(12|/V(ri,r2)|34)(-l)^+i-^(-l)j^+j3 

A:=0 A 

x5jJ'<5MM'((z4)jl||xf'^(l)ll(/3^)j3)((/2^)j2||x?'''^(2)||(Z4^)j4) 

(64) 

The 6-j symbols limit the possible allowed values of k and A. The reduced 
matrix elements are evaluated using [S7j and further limit k and A. Note 
that as the upper and lower Hartree spinors have different / values, the 
reduced matrix elements in Eqs. (|68j) and (|64j) must have the corresponding, 
appropriate I values. 

Now consider the remaining two-dimensional radial integrals, where the 
numbers are a shorthand for all the quantum numbers needed to uniquely 
specify the radial wave functions j54j . 

/■oo roc 

(12|/|(ri,r2)|34) = / / dridr2Ui(ri)U2(r2)/^(ri, r2)U3(ri)U4(r2) (65) 
Jo Jo 

Here R(r) = U(r)/r are the appropriate radial Dirac wave functions, in terms 
of Ga(r) and Fa(r), and again i = S,V. 

Using the Hartree spinor representation, the particle-hole matrix ele- 
ment is expressed as a sum of Dirac matrix elements of the types shown 
above [HSl 

<im = E(2J' + l){ ^j;:; j } (IbJ'lVlamJ') (66) 

No exchange term is required, due to the fact that the A and the nucleon 
are distinguishable particles here. For example, the particle-hole matrix 
element for the vector spatial component of the effective interaction is 

-^2;i4(vs) = (-l)J^+j^+^f:E(-l)'{ ^ ^ j]l /dridr2 

fc A 

X {Gi(ri)F3(ri)/feV(ri, r2)G2(r2)F4(r2) 

x((/lA^)jl||xf'^(l)ll(/3B^)j3)((/2A^)j2||xf'^(2)||(;4B^)j4) 
- Gi(ri)F3(ri)AV(ri,r2)F2(r2)G4(r2) 

x((/lA^)jl||xf'^(l)ll(/3B^)j3>((/2Bj)j2||xf'^(2)||(/4A^)j4> 



J h ji 
A js k 
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- Fi(ri)G3(ri)A^(ri,r2)G2(r2)F4(r2) 

x((/lB^)jl||xr^(l)ll(Z3A^)j3>((/2A^)j2||xr^(2)||a4B^)j4^ 
+ Fi(ri)G3(ri)/fe^(ri,r2)F2(r2)G4(r2) 

X ((/lB^)jl||xf'^(l)l|(^3A^)j3>((/2B^)j2||xf'^(2)||(/4A^)j4) 

(67) 

Here and are the I values corresponding to the upper and lower 
Hartree spinors respectively for the ith wave function where i = 1. . 4. Now 
the splitting, for a Si/2-doublet, is just the difference between the particle- 
hole matrix elements of the two available states, or 

X, _ „ J=jl+j2 _ J=|jl-j2| /fjnN 
— ^nA;nA ^nA;nA \^°) 

The substitutions used to acquire the appropriate indices for this case are 
n = 1, 3 and A = 2, 4. The solution to the Hartree equations yields a single- 
particle energy level for the GS, Ea- As previously mentioned, for the 
cases under consideration this level is in fact a doublet; however, Eq. (|^S)) 
evaluates only the size of the splitting. In order to determine the position 
of the doublet relative to E^, one needs the relation 

^ (2J + 1) 5e = (69) 
J 

We now have a framework with which to calculate the size of the Si/2" 
splittings of the single A-hypernuclei of interest here and to determine their 
location relative to Ea- The problem is reduced to Slater integrals and some 
algebra; the 6-j and 9-j symbols are determined using [621 163j . The Dirac 
wave functions needed to solve the radial integrals are taken as the solutions 
to the Hartree equations from the previous section. Once all the parameters 
in the underlying lagrangian are fixed, the splitting is completely determined 
in this approach as there are no additional constants fit to excited state 
properties jSH. We also mention that this approach is applicable to excited 
states and multiplets for this class of nuclei. 

To calibrate this approach, we apply it to ordinary nuclei. Two mod- 
ifications to our framework are required here. First, an exchange term is 
included because the proton and neutron are indistinguishable particles. As 
a result, the particle-hole matrix element becomes the following |53j 



^abilm 



i:(2j'+i){^» J } 

X [(lbJ'|V|amJ') - (-l)j^+j-+J'(lbJ'|V|maJ')] (70) 
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Experimental Data 


M2 Calculation 


GS E/B 


A ^ 


-11.69 ±0.12 




52 




-10.89 






-12.50 ±0.35 




47 




-12.03 




fSi 


-16.60 ±0.2 




49 




-17.37 




'S'2c 
A * 


-17.50 ±0.5 








-17.95 






-18.70 ± 1.1 




47 




-18.63 




^Pb 


-26.5 ±0.5 




49 




-27.81 


Eso 


A 


0.15 ±0.09 




50 




0.150 


Esp 


A ^ 


10.83 ±0.03 




50. 




8.849 






10.6 ±0.1 








8.314 






7.70 ± 1.0 




48. 




7.832 



Table 2. The experimental data used in the parameter fits. This in- 
cludes six GS binding energies (E/B), one spin-orbit splitting of the p- 
states (Eso = |Eipi/2 ~ Eip3/2 1)' three s-p shell A excitation energies 

(Esp — |Eip3^2 — Eisj^^l). The calculated values of these observables, using the 
M2 set, are also shown. These values are given in MeV. 



Second, the effective interaction is also modified, requiring additional cou- 
plings to the rho and pion fields [31] 



V(r,,) = 7lSf^ 



_„2 -msri2 ^ ^ . s 0-2 p-mvri2 

Isl _L ^(1)^(2) 

~ lU lU 



An ri2 



'^^ '^^ 4tt ri2 



(1) (2)i[^i:i[!^gi^:!^ 



11 IfM 



47r ri2 



^ ^ 47r ri2 



(71) 



These alterations make the ordinary nuclear matter case considerably more 
complicated than the case of single A-hypernuclei. 



5. Results 

5. 1 . Parameter Fits 

The full lagrangian contains a number of free parameters. Those con- 
stants which lie in the nucleon and meson sectors are fixed by the G2 pa- 
rameter set of FST given in Tabled In the A sector of the lagrangian, 
a total of six parameters remain undetermined to order = 3. Fits are 
conducted at various levels of truncation in the underlying A-lagrangian to 
fix the relevant constants. The fits performed here are entirely separate 
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Ml 


M2 


M3-1 


M3-2 


M4 


gSA/gS 


0.87357 


0.87697 


0.87390 


0.87090 


0.87697 


gVA/gV 


1.0 


0.98623 


0.97766 


0.98050 


0.98623 


gTA/4 




-0.892 


-0.891 


-0.877 


-0.892 


/^l 






-0.1550 


0.1500 


0.0700 


^2 






-0.2517 


0.2436 


0.3111 


^3 










0.0700 



Table 3. Lists of the constants in the five parameter sets constructed here. Note 
that all the constants are natural and that these sets represent different levels of 
sophistication in the A-lagrangian. 



from the one which determined the G2 parameter set; however, the frame- 
work which FST used to conduct their fits is identical to the one employed 
here. The experimental data utilized to constrain the parameters in the A- 
lagrangian is listed in Table [21 and consists of three types of observables: GS 
binding energies, s-p shell A excitation energies, and spin-orbit splittings of 
the p-states. Now we use the framework outlined in sections 21 and |31 to cal- 
culate these same observables for some initial guess of the parameters. The 
calculated and experimental values are both substituted into the equation 

i X 

where N is the number of data points and W^'' are the weights. The param- 
eters are varied such that the theoretical and experimental values converge. 
The constants are fixed at the values that produce a minimum in Xn- 

Our underlying A-lagrangian is truncated at four different levels and sep- 
arate parameter fits are conducted at each. First, we consider the simplest 
possible case; only terms to order v = 2 are retained in the A-lagrangian, 

(2) 

which corresponds to Cf^ ■ This A-lagrangian has a total of two free pa- 
rameters, gsA and gvA- In this case, the vector coupling is assumed to be 
universal, as it is coupled to the conserved baryon current, and the scalar 
coupling is fit to reproduce the binding energy of a single A in nuclear mat- 
ter, which is about -28 MeV |27]. These assumptions are in keeping with 
the previous work in The parameters determined here are shown in 

Table 131 as the Ml set. This set reproduces the GS binding energies fairly 
well, but is unable to simulate either the correct spin-orbit splitting in the 
p-states or the s-p shell excitation energies in light A-hypernuclei. 

In order to obtain a better fit to the data, we increase the level of 
truncation. Therefore, tensor couplings to both the vector and photon fields 



X, 



oxp 



X 



(72) 
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W 

c 

c 



i- (Ipyj) n(lP3/2)A 



— 3 -parameter Fit 

— Experimental Value 



10 



(ls,J, 



S -15 - 



-20 L 



0) 



10 - 



,C 



„(1Pi/2)a 



16 28 



3 -parameter Fit 

Experimental Value 



■20 - -1 



<D 

C " 

W 
to 
a 

•S -30 - 

m 



-40 L 



32 



40 



Ca 



208 



Pb 



A 



(NM) 



Fig. 1 . Results of the unweighted 3-parameter lit to a series of experimental data. 
The G2 parameter set of FST is used for both the nucleon and meson sectors Q]. 
The calculated binding energy of a single A in infinite nuclear matter is also shown. 





M2 


M3-1 


M3-2 


M4 


X^o(UW) X 100 


0.105 


0.0877 






X^o(W) X 10 


0.598 




0.515 


0.485 



Table 4. The values for both the unweighted (UW) and weighted (W) fits 
relative to the of the Ml set. Here is determined from Eq. (|72|l using 10 
pieces of data. 
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(T) 

are included, which correspond to the terms in C\ . As a result, a third free 
parameter, gxA, is introduced. This fit is performed using seven pieces of 
experimental data: the six GS binding energies and the spin-orbit splitting 
given in Table |5J In this particular case, the weights in Eq. H72|) are all 
taken to be equal. The resulting parameters are given in Table |21 as the M2 
set and all satisfy the assumption of naturalness. Table |21 also outlines the 
numerical results of this 3-parameter fit. The outcome of this fit is shown 
graphically in Fig.^ One can see that both the GS binding energies and the 
small spin-orbit splitting in the p-states are reproduced well. The calculated 
s-p shell excitation energies fail to duplicate the experimental values for the 
lightest A-hy per nuclei; however, it is correctly given by the time one gets to 
i°Ca. In Fig. HI the value of -32.4 MeV is given as the calculated binding 
energy of a single A in nuclear matter. This M2 parameter set will be used 
in the subsequent calculation of the Si/2-splittings. 

A plot of the proton, neutron, and A densities for the GS of ^'^Ca calcu- 
lated using this M2 set is shown in Fig. |21 A graph of the Hartree spinors 
from the A wave function, GA(r) and FA(r), for the GS of ^''Ca using the 
M2 set is also given in Fig. |21 Notice that the magnitude of the lower 
spinor is very small; this indicates that the A is essentially behaving as a 
nonrelativistic particle in the nuclear potential. 

Next, the two terms nonlinear in the scalar and vector field, shown in 
c}'J^\ are retained. This brings the number of unconstrained parameters up 
to five. For this 5-parameter fit, ten pieces of experimental data are used; 
in addition to the data utilized in the 3-parameter fit, the three s-p shell A 
excitation energies listed in Table I^J are also included. Two versions of the 
5-parameter fit were conducted here: one unweighted and one weighted. In 
the former case, all of the weights are equal. For the latter, the weighting 
scheme is as follows: W^' = 1.0 for GS binding energies; W^^ = 10.0 for 
s-p shell A excitation energies; and W]^ = 40.0 for the spin-orbit splitting. 

The weights were selected using the formula W^^ = /i(AEexp/Eexp) where 
fi is an arbitrary factor chosen to prevent any observable from dominating 
the fit [HSj. However, not enough similar data was available to constrain 
the two new parameters individually. As a result, we initially restrict these 
parameters with the relation 

^=f^V = 1.624 (73) 

where n.m. denotes the nuclear matter values [5]. The results of both 5- 
parameter fits are shown in Table El the M3-1 and M3-2 sets denote the 
unweighted and weighted schemes respectively. Again notice the parameters 
are all natural. However, the new parameters are not very well determined 
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and fail to significantly improve the fit in either case, as can be seen from 
Table m Therefore, we leave the constraint of Eq. (|73jl intact. 

Lastly, to include all possible terms in the A-lagrangian up to order 
v = 3, all three terms in C^^^ are retained. Again, not enough similar data 
was available to individually constrain the new parameters; therefore, we 
restrict these parameters with the relation 

^1 = ^3 = 0.225^2 (74) 

and fix the remaining constants using the M2 set. These ratios were chosen 
because they tend to concentrate the effects of the new contributions in 
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the surface of the nucleus, i.e. the additional contributions now vanish for 
uniform nuclear matter. This will have a greater effect on the s-p shell 
excitations than on the GSs. The weighting scheme described above was 
used. The resulting parameters are listed in Table |31 as the M4 set. Again, 
as seen in Table |l]the improvement in the overall fit is negligible. The M3-2 
and M4 sets both improve the fit to the GSs but do worse with respect 
to the s-p shell excitations; the M3-1 set has the opposite effect. Also we 
mention that the parameter sets M3-1, M3-2, and M4 yield very similar 
density distributions to those acquired from the M2 set. 

5.2. Sii2-splittings 

In this section we discuss the calculation of the Si/2-splittings in A- 
hypernuclei and the results obtained from these calculations. Following the 
methodology established in section |31 one needs to evaluate 5e from Eq. (|68|) 
to determine the size of these doublets. It is possible to separate (5e into 
contributions from each portion of the effective interaction, or 

8e = 5e{s) + (5e(vt) + 5e(vs) (75) 

where s, vt, and vs represent the scalar, vector time, and vector spatial com- 
ponents respectively. As it turns out, the scalar and vector time components 
each cancel in the splitting, shown by 

5€{s) = 5e(vt) = (76) 

Therefore, the Si/2-splittings are entirely determined from the vector spatial 
term in the effective interaction, or 

5e = (5e(vs) (77) 

This is true for any system in which either the A or the nucleon hole has 
j = 1/2. Note that this calculation tests a different sector of the underly- 
ing lagrangian than the mean field analysis and that, as there is no corre- 
sponding interpretation in the static limit (M — > cx)), it is here an entirely 
relativistic effect. Now, to determine the splitting we only need to evaluate 
the matrix element in Eq. (|67|) for the two appropriate J values. The in- 
tegrals are solved using the Hartree spinors, Ga(r) and Fa(r), calculated in 
the single-particle analysis. Notice that the integrals in the vector spatial 
contribution mix the upper and lower components of the Hartree wave func- 
tions. Numerically, the integration is performed using Simpson's method. 

The results of this analysis are contained in Table |S1 The splittings 
with a neutron hole listed in Table [5] all correspond to single-particle levels 
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Nucleus 


State 


Levels 


6e 




(lP3/2)pHlSl/2)A 


^GS' ^ 


426 




(1Pi/2)p(1si/2)a 




472 




(lP3/2)pniSl/2)A 


2i:l> 1- 


316 




(lpi/2)n(lSi/2)A 




480 




(lPl/2)n(lP3/2)A 




125 




(lPl/2)n(lPl/2)A 




661 




(ld5/2)r7^(lSl/2)A 


^GS' 2 


293 


A ^ 


(2si/2)n(lSl/2)A 


1gs> 0+ 


216 




(ld3/2)nHlSi/2)A 


2+ 1 + 
■^GS' ^ 


308 




(ld3/2)n^(lpi/2)A 




393 


TPb 


(lil3/2)n^(lSl/2)A 


7Ss> 6+ 


24 



Table 5. Si/2-splittings, and some excited states, are shown with their respective 
configurations, level orderings, and doublet magnitudes. Here LL denotes lower 
level and \Se\ is in keV. 



which were used in the fits of the preceding discussion, as shown in Fig. ^ 
The Si/2-splittings for }^0, j^Si, j^S, and ^'^Ca are plotted in Fig. |31 notice 
that these four splittings are all within the experimental error bars and 
that the appropriate level orderings are shown. It should be mentioned that 
the three excited states with neutron holes shown in Table 121 will overlap 
with other states of the same J value. Therefore in these cases one must 
diagonalize the hamiltonian to determine the correct splitting and level 
ordering. The remaining doublets in Table |SJ those with proton holes, are 
for predicted A single-particle levels. These three are shown in Fig. |2 here, 
in addition to the GS splittings for both and j^^N, the doublet for the 
first calculated excited state in N is also given. These splittings will be 
measured in an upcoming experiment using the reaction (e, e'K'^) with much 
greater resolution than the (7r"'",K"'") reactions |58l I59j. As the effective 
interaction used here is isoscalar, there is no distinction in this approach 
between proton and neutron holes. This is apparent when comparing the 
GSs of j^^N and ]^0; the slight difference in their splittings, which is only 
about 10 keV, arises from Coulomb effects. Also note that the splittings for 
configurations with the holes in the same shell are larger for the smaller j 
value. For example, the doublet for the GS of j^B, in the (lp3/2)p ^(1si/2)a 
configuration, is smaller than that of the GS of N, in the (lpi/2)p "'^(1si/2)a 
state. The level orderings for each calculated doublet are also given in Table 
13 Notice that for all of the cases considered here, the state with the higher 
J value is the GS or, in the case of excited states, the lower level. 
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Particle-hole Splitting 
Calculated Value 
Experimental Value 



M in 
iH -lU 

W -12- 
C -14 - 

-18- 
-20 
-22 ■ 



-2 - 



> 



£?-io 

W -12 
C -14 
C 



16 



-18 - 



16 



o 



28 



Si 



Particle-hole Splitting 
Calculated Value 
Experimental Value 



32 



40 



Ca 



Fig. 3. Graph of GS particle-hole splittings and their respective level orderings for 
]y^0 and ^*Si on the top and ^■^S and ^^Ca on the bottom. The single-particle 
calculations were conducted using the M2 parameter set and are plotted along- 
side the experimental values |^ 03 QHI ■ Notice that the splittings lie within the 
experimental error bars in all four cases. 



Recent gamma-ray spectroscopy experiments |57| (and the experimental 
error bars on the GS binding energy of j^^B) suggest that these particle- 
hole splittings are in fact much smaller. In addition, the measured GS 
spins of and j^^O are 1 and respectively [HniinHl! whereas the values 
predicted here is and 1 respectively. As the tensor coupling was important 
in the spin-orbit splittings, it may play an important role in the case of the 
Si/2-splittings. Higher order terms in the effective interaction, especially 
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Particle-hole Splitting 

-2 - — Calculated Value 
^ Experimental Value 

I "T (lP3/2)''p(l^./2)A ^. 



M -6 - — =~ 2 



c 

W .8- 

bo ,-1 
C /I- X 



c 
W 



-10 



-12- 



Particle-hole Splitting 
Calculated Value 



> Experimental Value 



2: 



c 1- !■ 



Experiment 



Fig. 4. Top: graph of particlc-holc splittings for \ B and \ N and their respective 
level orderings. In addition to the GSs, the first calculated excited state in J^^N 
is also included. These calculations were conducted using the M2 parameter set. 
The experimental value for the GS of j^^B is taken from [201 • Bottom: particle-hole 
splitting for the GS of fsPiy. The level orderings and splittings are shown for both 
theory and experiment. Here the G2 parameter set of FST was used pp. 



those involving the tensor coupling to the A, may be required to obtain a 
quantitative description of the small si/2-doublet sphtting and the correct 
level ordering. This is left for future work. 

The present analysis was also extended to the case of ordinary nuclei. 
The necessary modifications to the theory were discussed in section [l] We 
apply this approach to the case of 15P17 in the (2si/2)p(ld3/2)n state. As 
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noted before, this calculation will require direct and exchange contributions 
from the scalar, vector, rho, and pion terms in the effective interaction. 
Fortunately, the statement of Eq. (|76() holds here for the direct term and 
can be extended to include the direct rho time component as well. The 
result of our calculation is 413 keV; the observed value is 77 keV |6Uj . This 
is shown graphically in Fig. |5 notice that the correct magnitude and level 
ordering is obtained. However, it should be noted that this calculation is 
considerably more complicated than the A-N case. 

In summary, we have successfully extended the hadronic effective field 
theory developed by FST to the region of the strangeness sector correspond- 
ing to single A-hypernuclei. This framework has the intrinsic strength of 
directly incorporating the following: special relativity, quantum mechanics, 
the underlying symmetry structure of QCD, and the nonlinear realization 
of spontaneously broken chiral symmetry. Furthermore, DFT provides a 
theoretical justification for this approach. This lagrangian can be used 
for predictive purposes once all the free parameters are determined. As a 
result, it was of interest to make a minimalist extension of this methodol- 
ogy in which a single, isoscalar A is added to the theory. An appropriate 
A-lagrangian was constructed as an additional contribution to the full in- 
teracting lagrangian of FST. This system was solved using the Kohn-Sham 
analysis. Parameter fits were conducted at various levels of sophistication 
in the A-lagrangian while maintaining the full FST lagrangian with their 
G2 parameter set. The 3-parameter fit reproduces the GS binding energies 
and small spin-orbit splittings well, but fails to simulate fully the s-p shell 
excitations in the lightest hypernuclei, although by ^°Ca the correct exci- 
tation energy is obtained. The inclusion of additional parameters does not 
significantly improve the quality of the fit. 

Many of the GSs used in the fits were actually particle-hole states; as 
a result, it was of interest to calculate their splittings. A methodology for 
examining these splittings was developed using Dirac two-body matrix ele- 
ments of an effective interaction. This effective interaction followed directly 
from the underlying lagrangian and to lowest order corresponded to sim- 
ple scalar and vector exchange.El Note that this lagrangian was designed to 
calculate other phenomena and there is nothing contained in it that guar- 
antees the production of small particle-hole splittings. The primary con- 
clusion from the present analysis is that all of the results obtained for the 
Si/2"doublet splittings used in the fitting procedure lie within the current 
experimental error bars. As a partial calibration, a calculation of the GS 
particle-hole splitting in 15P17, a much more complicated case, achieved the 
correct level ordering and doublet magnitude. Using this approach predic- 
tions were made for nuclei that will be measured in an upcoming (e, e'K"^) 
experiment at Jefferson Lab |581 159j .l^ 
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In this appendix, we discuss the selection of the terms in our A-lagrangian 
to order v = 3. It is straightforward to see which terms are retained to order 
ly = 2, with the exception of the four fermion terms. Therefore, the following 
is a list of all remaining possible combinations of the fields to order v = 3, 
consistent with this approach, and a short discussion of each. 

• Four fermion terms in the nuclear case, such as NNNN, are eliminated 
by substituting the meson equations of motion into the lagrangian. 
Under normal circumstances this is not feasible; however, this is al- 
lowed when the system is already in equilibrium. Here we want to ex- 
tend the framework of FST to single A-hypernuclei with no additional 
mesons. In this case, either NNAA or AAAA can be eliminated using 
this method, but not both simultaneously. Fortunately, the second 
term involves self-fields of the A and consequently, can be discarded. 
This scheme also applies to terms with more than four fermion fields. 

• The term Ka^pW is consistent with this framework. 

• The terms KKcfP' and AAV^^ are consistent with this framework. In 
the nucleon sector, terms of this variety were regrouped using meson 
field redefinitions. Here the terms have different constants than in 
the nucleon case; therefore, these terms cannot simply be regrouped, 
unless additional mesons are included. 

• The term K'y^KcfN ^ is also retained. In the nuclear case, it was elim- 
inated via the Dirac equation, but this is not possible here. 

• Next, the following term is consistent with this methodology, but can 
be rewritten as 



Appendix A 




d 



(A.l) 



The second term is a total derivative, which does not change the ac- 
tion, and the third term is a four derivative of a conserved current, 
which is zero. Therefore this term can be neglected. 
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• Consider the following three terms: A7^0q|-A, A7^7jyV^g|-A, and 
A7ju7iyg|-^^A. The Dirac equation for the A can be substituted into 
each of these to convert them into a type of term already considered. 

• Lastly, all of the contributions with A^ are absorbed into other terms 
in the same manner as like terms with V^. However, the terms 
A7^AA^ and AAA^ can be discarded as Q = for the A. Therefore, 
the only remaining electromagnetic term is A(T^,yF^,yA. 

Note that the parameters have yet to be determined. When the terms are 
regrouped, the free parameters can be redefined to suit our purposes. 
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